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ABSTRACT 

We use a sample of 9397 low redshift (z < 0.1) galaxies with a close companion to investigate 
the connection between mergers and luminous infra-red (IR) galaxies (LIRGs). The pairs are 
selected from the Sloan Digital Sky Survey (SDSS) and have projected separations r p < 80 
kpc, relative velocities AV < 300 km s _1 and stellar mass ratios within a factor 1:10. 
A control sample consisting of 4 galaxies per pair galaxy is constructed by simultaneously 
matching in stellar mass, redshift and environment to galaxies with no close companion. The 
IR luminosities (L/«) of galaxies in the pair and control samples are determined from the 
SDSS - Infra-red Astronomical Satellite (IRAS) matched catalog of Hwang et al. (2010). Over 
the redshift range of our pairs sample, the IRAS matches are complete to LIRG luminosities 
(L/s> 10 11 L Q ), allowing us to investigate the connection between mergers and luminous 
IR galaxies. We find a trend for increasing LIRG fraction towards smaller pair separations, 
peaking at a factor of ~ 5-10 above the median control fraction at the smallest separations 
(r p < 20 IhjQ kpc), but remaining elevated by a factor ~ 2-3 even out to 80 ft™ kpc (the widest 
separations in our sample). LIRG pairs predominantly have high star formation rates (SFRs), 
high extinction and are found in relatively low density environments, relative to the full pairs 
sample. We also find that LIRGs are most likely to be found in high mass galaxies which have 
an approximately equal mass companion. We confirm the results of previous studies that both 
the active galactic nucleus (AGN) fraction and merger fraction increase strongly as a function 
of IR luminosity. About 7% of LIRGs are associated with major mergers, as defined within the 
criteria and mass completion of our sample. Finally, we quantify a SFR offset (ASFR) as the 
enhancement (or decrement) relative to star-forming galaxies of the same mass and redshift. 
We demonstrate that there is a clear connection between the ASFR and the classification of 
a galaxy as a LIRG that is mass dependent. Most of the LIRGs in our merger sample are 
relatively high mass galaxies (log (M*/M ) > 10.5), likely because the SFR enhancement 
required to produce LIRG luminosities is more modest than at low masses. The ASFR offers a 
redshift-independent metric for the identification of the galaxies with the most enhanced star 
forming rates that does not rely on fixed L/« boundaries. 
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1 INTRODUCTION 

The total far infra-red (IR) luminosity (L/ R ), which is denned as 
the bolometric contribution from 8-1000 yum, of a galaxy is closely 
tied to its star formation rate (SFR). The correlation is due to the 
processing of UV photons by galactic dust, which are re-emitted in 
the thermal spectrum (see, for example, the reviews on galactic star 
formation by Kennicutt 1998; Kennicutt & Evans 2012). There can 
be additional contributions to h 1R from an active galactic nucleus 
(AGN). However, most studies find that for the general population, 
the AGN contribution to the total h IR is modest, at about 10-20%, 
although the AGN contribution can become dominant at the highest 



far-IR luminosities (Nardini et al. 2010; Yuan, Kewley & Sanders 
2010; Petric et al. 201 1; Lee et al. 2012). 

Based on the two mechanisms for IR emission, it can be easily 
appreciated why mergers provide a fertile environment for high IR 
luminosities. It is well documented that mergers have, on average, 
enhanced SFRs relative to isolated counterparts (e.g. Kennicutt et 
al. 1987; Barton, Geller & Kenyon 2000; Bergvall et al. 2003; Lam- 
bas et al. 2003; Nikolic, Cullen & Alexander 2004; Alonso et al. 
2004, 2006; Li et al. 2008; Ellison et al. 2008, 2010; Woods et al. 
2010; Darg et al. 2010a; Patton et al. 2011). Although typical SFR 
enhancements are a factor of a few at both low and high redshift 
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(Robaina et al. 2009; Wong et al. 201 1), a minority of galaxies can 
have SFR enhancements of a factor of ten or more (Scudder et al. 
2012). Interactions can also lead to nuclear activity, triggering AGN 
over a wide variety of energies from Seyferts to quasars (Kennicutt 
et al. 1987; Ramos-Almeida et al. 201 1, 2012; Canalizo & Stockton 
2001; Ellison et al 201 lb; Silverman et al. 201 1; Koss et al. 2012). 
It is therefore not surprising that pairs also show increased L ;K and 
that interactions can contribute significantly to the IR luminosity 
(e.g. Lonsdale, Persson & Matthews 1984; Kennicutt et al 1987; 
Telesco, Wolstoncroft & Done 1988; Xu & Sulentic 1991; Xu et al. 
2000; Hernandez Toledo, Dultzin-Hacyan, & Sulentic 2001; Bridge 
et al. 2007; Smith et al. 2007; Hwang et al. 2010, 2011). 

At the highest IR luminosities, galaxies are classified as 
luminous IR galaxies (LIRGs) and ultra-luminous IR galaxies 
(ULIRGs) when their total far-IR luminosities are 1 1 < log 
(L ;R /L ) < 12 and log (L IR [L Q ) >12 respectively (see Sanders & 
Mirabel 1996 and Veilleux et al. 2002 for reviews). It is well known 
that, at least at low redshift, the vast majority of ULIRGs are merg- 
ers, with recent work suggesting that mergers may also contribute 
significantly to ULIRGs at high z (e.g. Dasyra et al. 2008; Kar- 
taltepe et al. 2010, 2012). As the h IR decreases, so does the frac- 
tional contribution from mergers (e.g. Kartaltepe et al. 2010), such 
that at log (L w /L ) = 11.5 about 20% of LIRGs are major merg- 
ers. Moreover, at fixed h IR , there is an increasingly important con- 
tribution from mergers at lower redshifts (Melbourne, Koo & Le 
Floc'h 2005; Wang et al. 2006). However, not all mergers lead to 
such prodigious IR emission. Statistically this is obvious from the 
comparison of the ULIRG space density of ~ 10~ 7 Mpc~ 3 at low 
redshift (Kim & Sanders 1998), compared to the space density of 
mergers (e.g. Chou, Bridge & Abraham 201 1) which is 3 orders of 
magnitude higher. 

In this paper, we study the IR properties of a large sample of 
low redshift close pairs in order to further investigate the connec- 
tion between high h IR and mergers. Whilst the fraction of LIRGs 
that are mergers has been previously quantified, the fraction of 
mergers that result in LIRGs is not well constrained. How fre- 
quently does a merger result in a LIRG, and can we identify which 
mergers are most conducive to LIRG production? We also compare 
the properties of the LIRGs in close galaxy pairs with LIRGs in 
our control sample, in order to understand whether merging plays a 
special role in LIRG production, or whether interactions are simply 
one pathway to prodigious IR emission. 

We adopt a cosmology of C1 A = 0.7, Q M = 0.3, H = 70 
km/s/Mpc. 



2 SAMPLE SELECTION 
2.1 Infra-red selection 

The study presented here combines a spectroscopically selected 
sample of pairs from the Sloan Digital Sky Survey (SDSS) with 
IR luminosities obtained from the Infra-red Astronomical Satellite 
(IRAS). IRAS surveyed the sky in four bands at 12, 25, 60 and 
100 yum. Hwang et al. (2010, hereafter H10) matched the IRAS 
Faint Source Catalogue (FSC), which contains 173, 044 IR sources 
at \b\ > 10 deg with the SDSS data release (DR) 7. H10 require 
that IRAS sources must have robust detections at 60 yum, but no re- 
quirements are placed on detections in the other bands. Bolometric 
IR luminosities are calculated by H10 via spectral energy distribu- 
tion (SED) fits to the templates of Chary & Elbaz (2001) at 60 and 
(where available) 100 yi/m. Hwang et al. (2010) show that IRAS 60 
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Figure 1. versus redshift for IRAS-detected galaxies matched to the 
SDSS DR7 in the Hwang et al. (2010) catalog. We adopt a redshift cut 
z <0. 1 (vertical dashed line) so that we are complete for LIRG luminosities 
of log (L/^/Lq) ^11 (horizontal dashed line). 



yum only determinations of h 1R agree well with those derived from 
60+100 yum SED fits, with a scatter of only 0.08 dex (see also Zhu 
et al. 2008). 

One of the main goals of this work is to study the properties 
and frequency of LIRGs in our galaxy pair sample. Figure [TJshows 
the distribution of redshift and L;« for galaxies in the H10 cata- 
log, and demonstrates the increasing incompleteness at small L w at 
larger distances. In order to be complete down to the limiting LIRG 
luminosity (log (L 1R /L e ) > 1 1), we impose a redshift cut z < 0.1. 

2.2 Close pairs sample 

We select a sample of close galaxy pairs from the subset of galax- 
ies in the SDSS DR7 whose extinction-corrected r-band Petrosian 
magnitudes are in the range 14.0 < m, < 17.77. For all galaxies, we 
compute local environmental densities, £„: 

where d„ is the distance in Mpc to the n' 1 ' nearest neighbour 
within +1000 km s -1 . Normalized densities, 6„, are computed rel- 
ative to the median Z„ within a redshift slice ± 0.01. We further 
restrict the parent sample to have a redshift range of 0.01 < z < 0.1. 
The lower redshift cut ensures cosmological values and the upper 
redshift cut ensures H10 catalog completeness for LIRGs. From 
this parent sample, close galaxy pairs are selected that fulfill the 
following criteria: 

(i) A projected physical separation r p ^ 80 kpc. 

(ii) A line of sight velocity difference AV <300 km s _1 , in order 
to minimize chance projections. 

(iii) A stellar mass ratio of 0.1 < M1/M2 < 10, as determined by 
Mendel et al. (2013fl 

(iv) A normalized local density estimator c5 5 < 75, i.e. the local 



In this paper, we consider the mass ratio to be a property of an individual 
galaxy, rather than the pair as a whole. That is, the mass ratio is the mass of 
a galaxy relative to its companion. This approach allows us to distinguish 
between the more (M1/M2 > 1) and less (M1/M2 < 1) massive galaxy in a 
pair. 
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density must be no more than 75 times more than the median den- 
sity in a redshift slice of ±0.01 around the galaxy. This criterion is 
imposed so that an adequate number of controls can be matched to 
the pairs sample, see below. 

The pair selection described above is very similar to that 
which we have used in previous papers (e.g. Patton et al. 2011), 
but with recently re-calculated (bulge+disk) masses from Mendel 
et al. (2013). As shown in Simard et al. (2011) and Patton et al. 
(2011), the SDSS pipeline photometry, which is used to deter- 
mine masses through spectral energy distribution (SED) fitting, 
can be seriously affected by the presence of a close companion. 
Mendel et al. (2013) have performed a thorough re-assessment of 
masses for the full spectroscopic SDSS DR7 using the improved 
(contamination-corrected) photometry of Simard et al. (2011). The 
star formation rates used in this paper are the 'total' (aperture- 
corrected) values determined from the SDSS spectra (Brinchmann 
et al. 2004; Salim et al. 2007). The star formation rates are de- 
termined through a Bayesian fitting procedure that uses spectral 
energy distribution models with a range of input parameters. For 
star-forming galaxies, the SFRs are essentially constrained using 
the strong Balmer emission lines, which agree well with UV mea- 
surements (Brinchmann et al. 2004; Salim et al. 2007). Aperture 
corrections are performed based on a colour dependent term of the 
light outside of the fibre region, as described in Brinchmann et al. 
(2004). 

The final step in the compilation of the pairs sample is a cor- 
rection for incompleteness at small angular separations due to fi- 
bre collisions. The relative excess of galaxies at wider separations 
(where the completeness rapidly increases) can be corrected by 
culling random pairs from the spectroscopic sample. This culling 
procedure randomly excludes 67.5% of pairs at angular separations 
6 > 55 arcseconds and has been shown in our previous work to 
effectively remove any bias in pair properties as a function of pro- 
jected separation (e.g. Ellison et al. 2008, 2011b; Scudder et al. 
2012). Of the 9397 galaxies in a pair which fulfill all of our selec- 
tion criteria and pass the cull, 400 are matched in the H10 catalog, 
which associates the closest SDSS counterpart to the IRAS detec- 
tion. 

Since the size of the IRAS beam is much larger than the 
galaxy-galaxy separations in our pairs, the h IR recorded in the H10 
catalog is the sum of the two galaxies. We discuss in Section [3~T| 
how we account for this in terms of classifying individual galaxies 
as LIRGs. More fundamentally, we can ask whether the galaxy- 
galaxy blending of IR flux is systematic for all the pairs in our 
sample, or whether some pairs will be resolved by the IRAS beam, 
introducing a separation dependent effect. The size of the IRAS 
beam is 1x5 arcmin, 1x5 arcmin, 2x5 arcmin and 4x5 arcmin at 
12, 25, 60 and 100 fim respectively. Most (95%) of our pair galaxies 
rely solely on the 60 j-im flux, so we will consider the effective IRAS 
beam size for the determination of L 1R to be 2x5 arcmin. 97.6% of 
our pairs sample have angular separations less than 2 arcminutes 
and 99.6% have separations less than 3 arcminutes. Therefore, for 
the vast majority of our pairs the total IR flux of the two galaxies 
is in the single beam. The L w measurements therefore consistently 
include the IR flux from both members of a paifl 

2 More recent IR satellites offer an improved spatial resolution, ranging 
from a few, to a few tens of arcseconds. However, without a careful treat- 
ment of deblending, the photometry of close pairs (r p < 20 Kji kpc) can 
be still severely affected by near neighbour contamination in imaging at 
arcsecond resolution (e.g. Simard et al. 201 1; Patton et al 201 1). It is there- 
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Figure 2. Distances to the 3rd and 5th nearest neighbour (top and bottom 
panels respectively) for pairs (red) and control (black) for galaxies that have 
been matched to the Hwang et al. (2010) IRAS catalog. For both density 
metrics there is excellent agreement between the pair and control distribu- 
tions. 



2.3 Control sample 

A matched control sample is an important component of under- 
standing the differential changes that occur in galaxies during the 
merger process (e.g. Perez et al. 2009). The sample of galaxies from 
which we construct our control is initially all galaxies which have 
no close spectroscopic companion within a projected separation of 
80 kpc and with relative velocity less than 10,000 km s _1 . How- 
ever, due to the spectroscopic incompleteness of the SDSS, this 
'control pool' will still contain some truly interacting galaxy pairs 
for which only one (or neither) of the galaxies has a measured red- 
shift. In order to remove contaminating mergers from the control 
pool, we make use of the visual classifications presented by the 
Galaxy Zoo project (Lintott et al. 2008). Darg et al. (2010b) de- 
scribe the identification of a merger sample from the Galazy Zoo, 
where they use an initial cut of the weighted merger vote fraction, 
f m > 0.4, followed-up by further visual confirmation. Although this 
vote fraction cut leads to a relatively pure merger sample, and al- 
lows a practical starting point for visual confirmation, Darg et al. 
(2010b) note that the first signs of visible disturbance can occur in 
galaxies with f„, as low as 0.2. Moreover, our own visual inspec- 



fore more straighforward to deal with the fully blended IRAS images, than 
with images that exhibit a 'sliding scale' of contamination as a function of 
projected separation. 
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Figure 3. L/« distribution for our control sample (black) and pairs (red). 
Galaxies with a close companion have systematically higher IRAS IR lumi- 
nosities. 



tion of low /„, galaxies reveals that many galaxies with f m ~ 0. 1 
may not be strongly disturbed, but can still have a close companion 
(indeed, we have close pairs in our spectroscopic sample with sim- 
ilarly low /,„), although the companion tends to be faint relative to 
the main galaxy. We hence conservatively exclude all galaxies with 
a non-zero f m from our control pool, which reduces its size from 
594,645 to 421,936 galaxies. 

The control sample used in this work is matched to the pairs 
sample simultaneously in mass, z and environment (local normal- 
ized galaxy density, 65). In previous work, we have taken an itera- 
tive approach of matching the best control galaxy without replace- 
ment until the control distributions are no longer consistent with the 
pairs (e.g. Scudder et al. 2012). This approach typically yields 10- 
12 matches per pair galaxy; a similar number is achieved for large 
samples that have their control samples only matched in mass and 
redshift (e.g. Patton et al. 2011). However, for a large pairs sam- 
ple, matching on three parameters yields only a small number (1 or 
2) of controls before the statistical comparison fails. The failure is 
usually in the distribution of densities; a small number of pairs with 
atypically high local densities have few possible matches amongst 
the controls. We therefore include a maximum density threshold in 
the pair selection criteria (Section 12.2) of 75 times the local me- 
dian. This excludes only ~ 2.5% of pairfl but allows us to match 4 
controls per pair galaxy. 

The environmental matching is particularly important for the 
present study, due to the large size of the IRAS beam. Even for our 
control galaxies there may be 'contamination' (additional IR flux) 
in the measured L;« from nearby galaxies. However, for a given 
environmental density, we expect the same contamination for the 
pairs and control. At z = 0.1 the largest dimension of the IRAS 
beam (5 arcmin) corresponds to more than 0.5 Mpc. Since pairs 
tend to be in moderately over-dense regions such as groups (Bar- 
ton et al. 2007), without environmental matching there may be ad- 
ditional IR flux from nearby galaxies (the contribution from the 
companion is dealt with separately, see Section [3~T| below). We can 
use our measurements of environmental density to gain some intu- 
ition over how well the density matching will perform on the scale 
of the IRAS beam. Figure [2] shows the distribution of distances to 



3 As we show in Section 13.2.21 most IRAS detected galaxies are in rela- 
tively low density environments. 



the 3rd and 5th nearest neighbours for the pair and control galax- 
ies that have been matched to the H10 catalog. Although one of 
the parameters in the control matching process was the locally nor- 
malized density, 65, we do not necessarily expect the pair and con- 
trol histograms in Figure |2]to be perfectly aligned. This is because 
only galaxies with IRAS counterparts are included, for both the 
pair and control distributions. That is, whilst the matching is done 
independently of IRAS detections, only those in the IRAS sample 
are shown in Figure|2] Nonetheless, in Fig.[2]fhere is good agree- 
ment between the distributions of ^5 in the pairs and control. There 
is also an excellent agreement between the distribution of the 3rd 
nearest neighbour distances of pairs and controls, as shown in the 
top panel of Figure|2](where, for a pair, its nearest neighbour is the 
companion). The environmental matching should therefore largely 
mitigate the issue of multiple galaxies within the IRAS beam. 

Figure [3] shows the L, 1R distribution of our pairs and control 
samples. The pairs are, on average, slightly more IR luminous than 
the full H10 sample, with a typical log (L 1R /L ) ~ 11. Indeed, of the 
400 paired galaxies which are matched with H10, 172 are classified 
as LIRGs (40%) and 1 as a ULIRG. It can already be seen from 
Figure [3] that the pairs have higher average L /R compared to their 
mass, redshift and environment matched controls. 



3 RESULTS 

3.1 The luminous infra-red galaxy fraction in close pairs 

We begin our investigation of the IR properties of pairs by quan- 
tifying in Figure [4] the fraction of LIRGs in our pairs sample as 
a function of projected separation. The LIRG fraction in each r p 
bin is calculated relative to the total number of pair galaxies in 
that r p range, regardless of IRAS detection. There is a steady in- 
crease in the LIRG fraction in the pairs as the projected separation 
decreases, indicating that their high h IR is associated with the in- 
teraction. In addition to the general increase in the LIRG fraction 
amongst the pairs sample, it is of interest to estimate the relative 
LIRG frequency compared to the control sample. The lower panel 
of Fig. [4] shows that the LIRG fraction in pairs is enhanced relative 
to the control sample at small separations (r p < 20 /i^J kpc) by a 
factor of ~ 5 - 10, and remains enhanced at the highest separations 
of the sample (80 hzl kpc) by up to a factor of ~ 2 - 3. 

It is important to consider whether the LIRG fraction enhance- 
ment is just extra IR flux moving into the beam as galaxies get 
closer together. Similarly, one might question whether the LIRG 
fraction plateau at wide separations is caused by residual environ- 
mental effects not accounted for in our control sample matching. 
As argued in Section [272] for all but the widest pairs (~2%), both 
galaxies are in the beam, so that the trend of increasing LIRG frac- 
tion as a function of r p can not be due to a gradual trend of beam 
confusion. Moreover, since the projected separations of this minor- 
ity of partially resolved galaxies are all r p > 30 kpc, separation 
dependent beam confusion/resolution can not explain the rise in the 
LIRG fraction at small projected separations. As a further test, we 
construct a sample of projected (in redshift space) pairs that we 
do not expect to show any signs of interaction, but are nonethe- 
less close on the sky. We have previously used such projected sam- 
ples to test for photometric contamination effects (e.g. Patton et 
al. 2011). The projected sample is selected identically to the main 
pairs sample, with the exception that the relative velocity criterion 
is now 300 < AV < 5000 kmr'.We first confirm that the SFRs of 
the projected sample do not depend on r p . Next, the LIRG fraction 
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Figure 4. Top panel: Fraction of galaxies in pairs (red) and control (open 
black circles) that are classified as LIRGs (L/« > 10" L Q ) as a function 
of projected separation. Lower panel: The ratio of LIRG fractions in pairs- 
to-control as a function of projected separation. The grey shaded circles 
have been corrected for the IR flux contribution from the companion galaxy. 
This figure reveals that the fraction of LIRGs is elevated in close pairs at all 
separations in our sample, averaging 5-10 times the control fraction at small 
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Figure 5. Pairwise properties for the full pairs sample (9397 galaxies, 
black), those detected by IRAS (400 galaxies, green) and those classified 
as LIRGs (172 galaxies, red). Mass ratios are log (M1/M2) i for the more 
massive companion and log (M1/M2) j for the lower mass companion in a 
pair. IRAS detections and LIRGs are biased towards more massive galaxies 
(see Figure|6] hence the asymmetry in the red and green histograms, relative 
to the full sample (black). 



of the projected pairs is calculated and found to be flat with r p , Had 
the IRAS beam been gradually missing flux at larger r p , then the 
LIRG fraction of the projected pairs would have decreased with in- 
creasing separation. We therefore conclude that the LIRG fraction 
of the AV < 300 km s pairs is truly declining with r p , and the 
effect shown in Fig.|4]is not a beam confusion effect. 

It should be noted that although our L ;R measurements are 
the combined fluxes from two galaxies, such blending is also fre- 
quently present in other LIRG and ULIRG samples (e.g. Howell et 
al. 2010). Nonetheless, it is of interest to attempt a decomposition 
of the hut between the two galaxy companions. As argued above, 
the environmental matching should largely account for this; sub- 
tracting the estimated IR flux from the companion is a very con- 
servative correction. We assume that the L ;R scales as the SFR, a 
relationship that is well-established in the literature and present 
in our control sample. The L ;s correction is therefore applied by 
dividing the total L; R between the two galaxies relative to their re- 
spective SFRs. Depending on the initial h IR in the H10 catalog, and 
the relative SFRs, this can either increase or decrease the number 
of LIRGs in each r p bin in Figure [4] For example, if a pair galaxy 
with log (L/g/Lg) = 11.1 that was originally classified as a LIRG 
has a companion with an identical SFR, the two galaxies will both 
be classified as non-LIRGs after the h IR has been equally divided. 



On the other hand, a pair galaxy with log (L/«/L Q ) = 11.9 would 
have originally been classified as a LIRG, but since no IRAS detec- 
tion is associated with its companion (due to the large IRAS beam), 
only 1 LIRG is counted for the pair. However, if the SFRs of the 
two galaxies are equal, the equal division of the L /R will yield two 
LIRGs. The overall result of the flux aportioning is shown by the 
grey shaded points in Figure |4] The SFR-corrected LIRG fractions 
generally decrease, but there is still a trend with r,, and an enhance- 
ment at all separations over the control by a factor of 2-10. 



3.2 Pairwise properties of LIRGs 

We next investigate whether there are certain pair configurations 
that are particularly conducive to producing LIRGs. Figure [^com- 
pares the properties of the pairs that are LIRGS, relative to the full 
pairs sample that fulfill the criteria listed in Section|2] The full pairs 
sample includes both actively star-forming and quiescent galax- 
ies. Since we expect the distribution of fundamental properties of 
these two categories to differ, we also compare the LIRG pairs with 
the pairs with IRAS detections, but with log (L /R /L ) < 11, which 
are therefore star-forming, but presumably at lower rates than the 
LIRGs, or with a lower dust content. 
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Figure 6. Galaxy masses and local galaxy densities (the environmental met- 
ric used in this work) for the full pairs sample (9397 galaxies, black), those 
detected by IRAS (400 galaxies, green) and those classified as LIRGs (172 
galaxies, red). LIRGs are preferentially found in lower density environ- 
ments relative to the full sample, but in a similar range of densities to other 
IR detected (star-forming) galaxies. 



3.2.1 Projected separation 

The top panel of Figure[5]shows that the LIRG pairs (red histogram) 
tend to have smaller projected separations than the full pairs sam- 
ple, as expected from Figure [4] A Kolmogorov-Smirnov (KS) test 
yields a probability of 1.6xl0~ 6 that the full pairs sample (black 
histogram) is drawn from the same r p distribution as the LIRGs (red 
histogram). This is also consistent with the finding by Hwang et al. 
(2010) that higher L /R galaxies tend to have closer nearest neigh- 
bours. The bias towards low r p is also seen in all the IRAS detected 
pairs (green histogram), where the LIRG and IRAS detected r p dis- 
tributions are indistinguishable in a KS test. This indicates that SFR 
enhancements contribute over a range of IR luminosities. There is 
also a hint of a deficit of LIRG pairs relative to the full pairs sample 
at wide separations, r p > 50 hzl kpc. Recall that the LIRG fraction 
relative to the control is still enhanced out to at least r p = 80 tCj^ kpc 
(Figure|4j. The relative paucity of LIRG pairs relative to non-LIRG 
pairs simply indicates that the SFR enhancements are likely to be 
more modest at these wide separations, as demonstrated explicitly 
in Scudder et al. (2012). 



3.2.2 Relative velocity 

The middle panel of Figure [5] shows that the LIRG pairs may be 
slightly skewed towards lower AV values than the full pairs sam- 



ple, although a KS test does not yield any significant difference 
between the AV distributions of any of the samples. Nonetheless, a 
bias towards smaller AV would be consistent with previous obser- 
vational results that have found that the SFR of close galaxy pairs 
increases with decreasing AV (e.g. Lambas et al. 2003; Nikolic et 
al. 2004; Alonso et al. 2004, 2006). However, the measurement of 
SFR dependence on AV is complicated by an underlying corre- 
lation between AV and local galaxy density (Ellison et al. 2010), 
which means that AV is itself a coarse metric of environment. In 
turn, there is also a well established SFR-density relation in the lo- 
cal universe (e.g. Gomez et al. 2003), so that a bias towards small 
AV values in LIRGs could be associated with higher intrinsic SFRs 
in lower density environments. In the bottom panel of Figure[6]we 
show that the LIRG pairs are certainly biased towards lower den- 
sity environments than the full pairs sample. However, the density 
distribution of LIRGs is similar to that of other IRAS detected star- 
forming pairs (green histogram in Figure [6] (see also Hwang et 
al. 2010). The environments of the full pairs sample approximates 
a double log-normal density distribution, where the high density 
population is apparently dominated by quiescent galaxies (Baldry 
et al. 2006). 

3.2.3 Mass ratio 

Our complete pairs sample is, by design, symmetric in mass ratio 
(black histogram in Figure |5}- However, the pairs that are identi- 
fied as LIRGs are biased towards high mass ratios (red histogram 
in Figure|5}, i-e. the LIRG galaxies in pairs tend to be more massive 
than their companions. A KS test yields a 2.3xl0~ 5 probability that 
the LIRGs are drawn from the same distribution of mass ratios as 
the full pairs sample. There are very few LIRG galaxies that have 
masses less than a factor of 3 below that of their companion. The 
magnitude limited nature of the SDSS survey means that high mass 
ratio galaxies are amongst the most massive in the sample. The ten- 
dency for LIRG pairs to have high mass ratios is therefore likely to 
be a selection effect. Galaxies with higher masses (and intrinsically 
higher SFRs) require a smaller SFR boost to achieve LIRG status 
(a topic which is discussed in more detail in Section |4~3V Indeed, 
relatively few of the LIRGs in either our pair or control samples 
have stellar masses log (M t /Mg) < 10 (Figure [6}. However, the 
lower panel of Figure [5] hints that the LIRG pairs have a narrower 
distribution of mass ratio than the IRAS detected pairs (although a 
KS test is inconclusive), a result which could not be explained by 
mass bias. Simulations have previously predicted that the most ef- 
fective SFR boost occurs within almost equal mass mergers (Cox et 
al. 2008). Observational works that have measured direct enhance- 
ments in the SFRs of close pairs support the dominant role of major 
mergers in producing the strongest starbursts (Woods et al. 2006; 
Ellison et al. 2008; Scudder et al. 2012; Lambas et al. 2012). The 
largest IR luminosities may therefore be associated with galaxies 
that are both massive (and therefore have relatively high SFRs to 
start with) and are engaged in major mergers (Hwang et al. 2010). 
Indeed, massive galaxies engaged in major mergers seem to domi- 
nate the local ULIRG population (Dasyra et al. 2006). 

To further investigate the dependence of LIRG classification 
on mass and mass ratio, in Figure [7] we plot the median h /R as a 
function of the total stellar mass of the pair. We distinguish the 
galaxy pairs by their relative masses, in order to test whether, at 
fixed total mass, the mass ratio has any impact on the total L/« (i.e. 
SFR). Galaxies with masses within a factor of three of each other 
are considered to be equivalent to major mergers (shown in red in 
Figure [7}; all others (up to a maximum mass contrast of a factor 
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Figure 7. h/ R as a function of the total pair mass (M[ +M2). Red points indi- 
cate major pairs (0.33 < Mi /M2 < 3) and black points show minor pairs up 
to the maximum mass ratio of our sample, 1:10. The horizontal dashed line 
is drawn at log (L/k/Lq) = 11. The offset to higher L/k at fixed total stellar 
mass in major pairs indicates higher levels of triggered star formation. 



of 10) are designated as minor mergers (black points). Figure [JJ 
shows that h /R increases with the total stellar mass of the pair, as 
expected from the correlations of mass with SFR and h IR . However, 
it can also be seen that at fixed total stellar mass, the median h IR 
is systematically higher for major mergers. The systematic offset 
to higher h IR at fixed total stellar mass in major mergers supports 
higher rates of dusty star formation in approximately equal mass 
encounters. 



3.3 The role of AGN 

In the previous sections, we have tacitly assumed that the main 
driver of the LIRG phase in our pairs sample is the enhancement 
in the SFRs that is a well-established result of galaxy interactions. 
In moderate luminosity LIRGs, the h IR indeed appears to be dom- 
inated by star formation (e.g. Alonso-Herrero et al. 2010). How- 
ever, AGN can also contribute to the IR luminosity. Petric et al. 
(201 1) confirm the modest contribution from AGN in the majority 
of LIRGs from a study of high ionization lines in the mid-IR spec- 
tra of 248 nearby LIRGs. They find that ~ 20% of LIRGs contain 
an AGN, and that the AGN is the dominant contributor of IR flux in 
only 10% of LIRGs. The relative importance of AGN to the IR flux 
increases with L ;R , so that the ULIRGs exhibit both high AGN frac- 
tions, and higher contribution from AGN within a given galaxy (e.g. 
Yuan et al. 2010 and references therein). In isolated star-forming 
galaxies, there is a well-known increase in the AGN (as classified 
from optical emission line ratios) fraction with stellar mass. Since 
the SFR also increases with stellar mass, the connection between 
AGN incidence and h IR is not too surprising. However, there seems 
to be an additional connection between the AGN incidence in IR lu- 
minous galaxies and mergers. The AGN-dominated LIRGs (Petric 
et al. 2011) and ULIRGs (Yuan et al. 2010) seem to preferentially 
occur at late merger stages. Combined with the observation that 
our pairs show an enhanced AGN fraction towards small separa- 
tions (Ellison et al. 2011b), it is germane to consider whether the 
increase in LIRG fraction seen in Figure|4]is due to an increasingly 
dominant contribution from AGN. 

Figure [8] shows the h IR distributions for pairs and control (red 
and black histograms, respectively) in star-forming, composite and 
AGN galaxies. The three categories are defined by combining the 
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Figure 8. h/ R distributions for star-forming (top panel), composite (middle 
panel) and AGN (lower panel) galaxies in our close pairs (r p < 30 kpc, 
red) and control (black) samples. The enhancement of in pairs relative 
to control is seen in the star-forming and composite galaxies, but not in the 
AGN. 



classifications of Kewley et al. (2001) and Stasinska et al (2006), 
and are made for galaxies with emission line signal-to-noise ratios 
S/N>3. For this figure, we consider only pairs (and their controls) 
with r p < 30 kpc, where differences in SFR, AGN and LIRG 
fraction (relative to the control sample) are most significant. Recall 
that we are only complete for log (L/«/L ) >11, but this incom- 
pleteness applies equally to all our galaxy classes. It is also im- 
portant to keep in mind that amongst the emission line galaxies in 
our sample, there is a higher AGN fraction at higher stellar masses, 
where the SFRs (and hence L IR ) are, on average, higher. Figure [8] 
shows a shift to higher h 1R for AGN versus star forming galaxies, 
e.g. by comparing the control samples for the two categories (black 
histograms in lower and upper panels respectively). Such a shift 
is expected based on the mass distributions of the AGN and star- 
forming galaxies, as described above. Figure |8]also shows that the 
pairs (red) have systematically higher h IR than the mass-matched 
controls, for both star-forming and composite classifications. The 
mass-matching means that the L IR differences between the pair and 
control galaxies is driven by the interaction, rather than a mass bias. 
However, there is no statistically significant difference between the 
L/r distributions of the pair and control AGN samples (Figure [8] 
lower panel). A KS test yields a ~ 50% probability (i.e. inconclu- 
sive) that the two samples are drawn from different populations. We 
conclude that the increase in LIRG fraction seen in Figure |4]is not 
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Figure 9. Fraction of galaxies in close pairs (r p < 30 /ujj kpc, red) and 
control (open circles) that are classified as AGN as a function of hi R . 

due to an increasing contribution to the h IR from AGN, but rather 
is driven by star formation. 

Having argued that the LIRG pairs are dominated by star for- 
mation, we can now turn to the question of whether the merger con- 
tributes to the hut (or mass) dependence of AGN incidence. Figure 
[9]shows the fraction of AGN as a function of L ;fi for both our close 
(r p < 30 \C.l kpc) pairs and control sample. Again, we recall that 
although our sample is only complete to log (L /R /L G ) > 11, the 
pairs and control can be compared differentially. In agreement with 
previous studies (e.g. Goto 2005; Hwang et al. 2010; Kartaltepe et 
al. 2010; Yuan et al. 2010; Iwasawa et al. 201 1) the AGN fraction 
increases steadily in the control with increasing h IR . As described 
above, this is likely due to the co-dependence of stellar mass on 
both SFR (which correlates with L /R ) and AGN fraction, such that 
higher stellar mass galaxies have higher SFRs (and hence higher 
L/ R ) and higher AGN fractions. Note that the exact AGN fraction 
at fixed h IR is highly dependent on the classification process (e.g. 
different line ratio cuts in the optical, colour classifications in the 
mid-IR and X-ray classification), making a direct comparison with 
other works non-trivial. Even for a given technique, such as optical 
line ratio diagnostics, precise classification is sensitive to choices 
such as S/N, as lower S/N cuts will increasingly include low ion- 
ization nuclear emission line regions (LINERs, e.g. Cid-Fernandes 
et al. 2010). The AGN classification of Stasinska et al. (2006) in- 
cludes composites in which even a minor contribution of AGN is 
inferred from the line ratios. Using a 'stricter' AGN classification, 
such as that of Kauffmann et al. (2003) reduces the AGN fractions 
by a factor of about 2, but still leaves the trend with h IR intact, 
yielding an AGN fraction of ~ 20% for 10 < log (L, R /L e ) < 1 1 and 
about 50% for the LIRGs. 

Beyond the general increase of AGN fraction with L IR , Figure 
[9]hints that the AGN fraction in pairs may be slightly lower than the 
control at fixed IR luminosity. Such behaviour might be expected 
if pairs increased their L IR without increasing their AGN fraction. 
This is in apparent disagreement with several recent studies that 
have linked mergers with AGN triggering in the local universe (El- 
lison et al. 2011b; Liu et al. 2012; Koss et al. 2012). One solution 
to this tension may be that the dust that is producing the high IR 
luminosities in our sample is hindering the classification of optical 
AGN. A forth-coming paper which uses multi-wavelength mid-IR 
diagnostics will investigate dust-obscured AGN in close pairs. 



4 DISCUSSION 
4.1 Merger fractions 

It is widely accepted that the fraction of mergers increases as a 
function of L /R , with mergers representing the bulk of the ULIRG 
population (e.g. Kartaltepe et al. 2010). At LIRG luminosities, 15- 
50 % of galaxies are usually classified as mergers (e.g. Zheng et 
al. 2004; Melbourne et al. 2005; Wang et al. 2006; Kaviraj 2009; 
Kartaltepe et al. 2010). However, these studies all start with an IR 
selected sample and classify mergers based on visual morpholo- 
gies. Our study of spectroscopically selected close pairs can offer 
an alternative measure of the merger fraction as a function of h IR . 

In order to determine the fraction of mergers as a function of 
L/jj, we must first define the term 'merger'. As for the AGN statis- 
tics in the previous section, we set an upper r p limit of 30 h^l kpc 
as the separation within which the strongest interaction-induced ef- 
fects are seen. The parent pairs sample has a mass ratio range from 
1:10 to 10:1. However, at higher redshifts the magnitude limited 
nature of the SDSS means that minor mergers become increasingly 
scarce and we are dominated by pairs with approximately equal 
masses (and magnitudes), see Section [3.2.3l It is additionally nec- 
essary to consider the mass completeness of our sample. For galax- 
ies near the (mass) detection threshold of the SDSS, we will be 
unable to detect lower mass companions. 

In order to tackle the above issues of completeness, the mea- 
surement of a merger fraction requires a careful definition of 'close 
companion'. The first component of this definition is in terms of 
companion proximity, both in projected and radial space. We adopt 
the criteria of a close companion to be a galaxy within 30 fci kpc 
and 300 km s -1 , We also restrict our definition to major mergers, 
i.e. galaxies whose companion has a mass within a factor of three 
of its own. The mass completeness of the SDSS is strongly colour 
dependent . At fixed redshift, blue galaxies can be detected to con- 
siderably lower masses than red galaxies (e.g. van den Bosch et al. 
2008; Mendel et al. 2013), due to the flux limit in the SDSS. At 
a redshift of z = 0.1, the SDSS is complete for both red and blue 
galaxies above log (M*/M Q ) ~ 10.7 (Mendel et al. 2013). A com- 
plete census of all 0.33 < Mi/M 2 < 3 companions would therefore 
require a mass limit at z = 0.1 of log (M*/M ) > 11.2, and hence 
severely limit the sample for which we could determine merger 
fractions. We therefore adopt a mass limit that only includes blue 
cloud galaxies (Mendel et al. 2013), and caution that our close pair 
fractions are only applicable in this mass range (at z = 0. 1, the mass 
limit corresponds to (M+/M ) ~ 10.37; for lower redshifts, the limit 
will be lower). 

The process for the quantification of merger fraction as a func- 
tion of L ;s is as follows. For each galaxy in the SDSS spectroscopic 
sample that appears in the Hwang et al. (2010) IRAS catalog, we 
determine whether its stellar mass is at least 0.5 dex (a factor of 
three) above the redshift dependent completeness for blue galax- 
ies described in Mendel et al. (2013). Every IRAS-detected SDSS 
galaxy that passes this mass criterion is then compared with the 
subset of our catalog of close pairs that have r p < 30 hzi kpc, AV < 
300 km s' 1 , and 0.33 < M,/M 2 < 3 to identify whether the galaxy 
has a companion. During the compilation of the pairs sample, pairs 
with angular separations 9 > 55 arcseconds are randomly culled 
to produce an even fibre (in)completeness as a function of separa- 
tion. Patton & Atfield (2008) showed that the level of spectroscopic 
completeness at these small separations is ~ 26%, relative to the 
photometric galaxy catalog. This completeness rate includes both 
the effects of fibre collisions and the overall spectroscopic com- 
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Figure 10. Fraction of galaxies that have a close companion within 30 ft^J 
kpc, 300 km s _I and 0.33 < M1/M9 < 3 as a function of L/r. All merger 
fractions have been corrected for fibre collision incompleteness and the gen- 
eral spectroscopic incompleteness of the SDSS. 
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pleteness of the SDSS. We therefore multiply our pair fractions by 
a factor of 3.85. 

Figure [Tol shows the fraction of galaxies in each h IR bin that 
have a close companion within the criteria for proximity and mass 
completeness described above. We find that below log (L; R /L G ) ~ 
1 1 the fraction of galaxies that have a close, approximately equal 
mass companion is approximately constant at ~ 2%. Above LIRG 
luminosities there is a rise in the fraction of galaxies that have a 
close companion, with an overall pair fraction amongst the LIRGs 
of 7%. It is not clear why the pair fraction increases only above log 
(L ;fi /L ) > 11, but a similar dependence has been previously seen 
(e.g. Kartaltepe et al. 2010). Tekola, Vaisanen & Berlind (2012) 
have shown that for galaxies with log (L IR fL Q ) < 11 there is no 
correlation between IR luminosity and local galactic density. How- 
ever, they also showed that for LIRGs, there is a strong correlation 
between h IR and density, and suggested that the LIRG luminosity 
(at low redshifts) represents a real transition in galaxy properties. 

Although the close pair fraction at LIRG luminosities that we 
find is of a similar magnitude to what has been found in previous 
studies (e.g. Zheng et al. 2004; Melbourne et al. 2005; Wang et al. 
2006; Kaviraj 2009; Kartaltepe et al., 2010), a rigourous quantita- 
tive comparison (i.e. on the basis of absolute merger fractions) with 
other works is non-trivial. For example, we have defined our merg- 
ers in a specific projected separation, AV and mass ratio range. The 
fraction of LIRGs (and SFRs of pairs in general) remain elevated 
beyond our nominal cut-off of r p < 30 hzl kpc. The mass ratio 
cut is dictated by completeness considerations, but minor mergers 
can experience significant SFR enhancements (e.g. Scudder et al. 
2012) that can result in LIRGs (Figure [5}. Indeed, Hwang et al. 
(2010) find that -34% of the LIRGs in their low redshift sample 
(also a result of a cross-matching between IRAS and SDSS) are in 
minor mergers. Our close pair fraction of 7% might therefore be 
considered to be a lower limit to the actual fraction of LIRGs with 
a companion. Conversely, our close pairs sample will also include 
galaxies that have not yet interacted (i.e. have not yet experienced 
a pericentric passage), as well as some projected pairs that are not 
actually destined to interact. Both of these categories will dilute the 
fraction of LIRGs associated with merger-triggered star formation. 
Scudder et al. (2012) visually classified a subset of our star-forming 
pairs and found only ~ 55% to show visible signs of disturbance. 
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Figure 11. Top panel: The distribution of SFRs in control galaxies for 
LIRGs and non-LIRGs. Middle panel: The distribution of SFRs in pair 
galaxies for LIRGs and non-LIRGs. Bottom panel: The distribution of SFRs 
in LIRG galaxies in the pair sample and control sample. In all panels, pairs 
are shown in red, control in black; LIRGs are shown with a solid line and 
non-LIRGs with a dot-dashed line. 

Finally, we do not include a merger timescale, or an observability 
timescale. These various caveats serve to caution about an abso- 
lute comparison with other works, but the general trend of a rising 
close pair fraction at high h /R should be robust, and the results in 
Figure [TO] demonstrate a connection between mergers and high IR 
luminosities. 

4.2 Star formation rates and extinction 

In Section [331 we argued that the increase in LIRG fraction was not 
dominated by the increasing AGN fraction in our close pairs sam- 
ple (Ellison et al. 201 lb). The enhancement in the SFR in pairs is 
therefore the obvious culprit responsible for the increasing LIRG 
fraction at small separations. In Figure[TT]we show the distribution 
of SFRs for the pair and control samples, distinguishing between 
those that are classified as LIRGs and those with lower h IR . In the 
top and middle panels we compare the LIRG (solid lines) and non- 
LIRG (dot-dashed lines) SFRs for control (top) and pairs (middle). 
It can be seen that the SFRs of the LIRGs are shifted to higher 
values, in both the pairs and control, relative to the SFRs in the 
non-LIRGs. In the lower panel, we compare the SFRs of the pair 
LIRGs (red) with the control LIRGs (black). The SFRs are indistin- 
guishable; the controls that are LIRGs have similarly high SFRs to 
those pairs that are LIRGs. Figure[TT]indicates that whilst the inter- 
action serves to increase the SFR, and hence the L 1R , similarly high 
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Figure 12. Top panel: The distribution of E(B - V) in control galaxies for 
LIRGs and non-LIRGs. Middle panel: The distribution of E(B - V) in pair 
galaxies for LIRGs and non-LIRGs. Bottom panel: The distribution of E(B— 
V) in LIRG galaxies in the pair sample and control sample. In all panels, 
pairs are shown in red, control in black; LIRGs are shown with a solid line 
and non-LIRGs with a dot-dashed line. 



values of SFR (and L;«) can be found amongst the control sample. 
This result is in contrast with the ULIRGs, where the high SFRs 
are almost exclusively associated with mergers. 

Star formation alone does not produce high IR fluxes; it is the 
presence of dust to re-process the shorter wavelengths that is the 
second key ingredient in the production of a LIRG. In Figure[T2lwe 
present distributions of E(B - V) for the LIRG/non-LIRG pairs and 
control in the same fashion as for SFRs in FigureQT] The E(B - V) 
values are determined from the Balmer decrement using an SMC 
extinction curve (see Scudder, Ellison & Mendel 2012 for full de- 
tails of our emission line processing). As was the case for SFRs, 
we find that the LIRGs in both the pair and control samples have 
higher E(B-V) than the non-LIRGs. Moreover, the pair LIRGs and 
control LIRGs (lower panel) have similar E(B - V) distributions. A 
high extinction may be due to an intrinsically higher dust content, 
or it can be an inclination effect. We investigated the disk incli- 
nations derived for our sample from the public SDSS photometric 
catalogue of Simard et al. (201 1), but found no difference between 
the LIRG and non-LIRG distributions. The high ~L IR in LIRGs is 
therefore apparently not an inclination effect, but relies on intrinisi- 
cally high SFRs in the presence of dust. 



4.3 How to make a LIRG 

It is clear from the results presented in this paper, and elsewhere 
(e.g. Robaina et al. 2009; Kartaltepe et al. 2010; Hwang et al. 2010, 
2011), that mergers contribute to the LIRG population. In this pa- 
per, we have demonstrated the LIRG-merger connection in two 
main ways. First, we have shown that the LIRG fraction increases 
in our sample of close pairs as the galaxies' projected separation 
decreases (Figure|4j. Second, we have shown in Figure [10] that the 
contribution from mergers increases at IR luminosities above the 
LIRG cut-off. We have also demonstrated that equal mass pairings 
of massive, star-forming galaxies are particularly conducive to pro- 
ducing LIRGs. However, as demonstrated in Figure [10] mergers 
are clearly not necessary for the production of LIRGs; any process 
which enables the galaxy to reach a high enough SFR (in the pres- 
ence of dust) may lead to a LIRG. This logic implies that there is 
some SFR threshold required for a galaxy to appear as a LIRG. In- 
deed, as shown in FigurefTTIthere are very few LIRGs with SFR < 1 
M Q /yr. This is not surprising, given the strong correlation between 
h IR and SFR (e.g. Kennicutt 1998; Kennicutt & Evans 2012). Form- 
ing a LIRG could therefore be as simple as achieving sufficient SFR 
(in the presence of dust). 

We can approximate^ determine the SFR threshold required 
to produce log (L 1R /L e ) > 11 empirically from our data by per- 
forming a linear least squares fit between the L, 1R values from the 
H10 catalog and the SFRs from the matched SDSS spectroscopic 
sample. We find that log (L /R /L ) =11 corresponds to a SFR ~ 9 
M /yr. Alternatively, from the L /R -mass relation (that results from 
correlations between SFR-L /;; and SFR-mass) we find that a log 
(L/k/Lq) =11 corresponds to a typical mass log (JVL/M ) ~ 10.9, 
at z ~ 0. 'Normal' star-forming galaxies above log (M*/M ) = 
10.9 may therefore readily manifest themselves as LIRGs, modulo 
the observed scatter in the stellar mass-SFR relation for individual 
galaxies. Lower mass galaxies require a SFR enhancement above 
the average for their stellar mass. This enhancement is simply the 
difference between a SFR of 9 M Q /yr and the value predicted from 
the stellar mass-SFR relation and is determined to be 



ASFR = 8.709 - 0.797 log M*. 



(2) 



The quantity ASFR is identical to the SFR offsets that we have 
previously calculated for the galaxy pairs in Scudder et al. (2012). 
That is, for a given galaxy in a pair, the difference between its SFR 
and the median of its mass-matched controls. For the present work, 
we have expanded the calculation of ASFRs beyond just the pairs 
sample, to include all galaxies with measured masses in the DR7 
spectroscopic sample (i.e. the sample presented in Mendel et al. 
2013). The ASFR of a given galaxy is calculated relative to all 
galaxies on the star-forming 'main sequence' (i.e. the SFR-stellar 
mass relation) within a redshift tolerance Az=0.005 and mass tol- 
erance Alog M*=0.1 dex. At least 5 matches are required in or- 
der to calculate a SFR offset, a limit which is achieved in 93% of 
cases (the typical number of mass and redshift matches is usually 
of the order several hundred). If less than 5 matches are found, 
the tolerance is grown by a further Az=0.005 in redshift and Alog 
M*=0.1 dex in stellar mass until the required number of matches is 
achieved. 

With SFR offsets in hand, it is possible to directly compare the 
ASFR of the LIRG pairs with the relation in equation [2] In the top 



The finer details of a galaxy's L;« will depend on factors other than just 
a SFR. Dust must obviously be present, and the spatial distribution of stars 
and dust may also play a role. 
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panel of Figure Q~3] we plot the ASFRs of the pairs vs their stellar 
mass; large circles indicate galaxies identified as LIRGs (after the 
flux is aportioned relative to SFR) and small points show non-LIRG 
galaxies. The diagonal dashed line shows the minimum ASFR re- 
quired in order to produce a LIRG (from equation |2]. The lower 
panel also depicts the ASFR versus mass, but for the control galax- 
ies. Despite its simplicity (and the significant scatter in the h IR - 
SFR relation, and smaller scatter in the stellar mass-SFR relation) 
equation [2] is a fairly good representation of the data, with most of 
the LIRGs lying above the diagonal dashed line in both panels of 
Fig. HU and most non-LIRGs lying at smaller values of ASFR for 
their stellar mass. Nonetheless, there are some non-LIRG galaxies 
whose A SFR appears, at face value, to be sufficient to produce a 
LIRG. The non-LIRGs in the shaded regions of Figure[T3]do have 
E(B-V) values that are lower, on average, than the LIRGs, but there 
is significant overlap. Presumably, the location of the dust relative 
to the starburst is important for the production of large IR fluxes. 

It can be seen that, in both the pairs and control samples, 
LIRGs are present across the full mass range, but larger ASFRs 
are required at lower masses, by a factor of 10 or more, as pre- 
dicted from equation|2] SFR offsets of a factor of 10 are very rare; 
Scudder et al. (2012) find only 5% of pairs have log A SFR > 
1. Figure [13] shows that SFR offsets of a factor of 10 are com- 
pletely absent from the control sample. The fact that LIRGs are 
more common at higher mass (e.g. Melbourne et al. 2008; Hwang 
et al. 2010), is presumably due, at least in part, to the lower ASFR 
that is required to put them above the LIRG threshold (diagonal 
dashed line in Figure [13}. Put another way, it is not that mergers 
between massive galaxies necessarily produce the largest SFR en- 
hancements (in fact, Figure [13] shows a spread of ASFR at fixed 
mass), but rather they require relatively small boosts (or no boost at 
all!) in order to achieve a LIRG-sufficient SFR. Indeed, the fraction 
of LIRGs with peculiar morphologies (presumably associated with 
an interaction) increases towards lower redshifts (Melbourne et al. 
2005). Secular processes, such as bars, which can enhance SFRs 
by similar amounts to mergers (e.g. Ellison et al. 2011a and refer- 
ences therein) may play a significant role in producing low redshift 
LIRGs (e.g. Wang et al. 2006). 

4.4 ULIRGs 

In this paper, we have focused on LIRGs, largely because they 
are the highest h IR galaxies we can study with reasonable statis- 
tics in our sample. The one ULIRG in our sample is IRAS 
F08572+3915 (SDSS objID=588013382191349914). This is a 
well known ULIRG and is the only galaxy pair in our sample that is 
present in the IRAS Revised Bright Galaxy Sample (Sanders et al. 
2003). This galaxy pair is also present in the Great Observatory All- 
sky LIRG Survey (GOALS, Armus et al. 2009). The IR luminosity 
for this ULIRG (from Hwang et al. 2010) is log (L«/L ) = 12.01, 
slightly lower than reported in the catalog of Wang & Rowan- 
Robinson (2009) log (L /fl /L G ) = 12.16. If we aportion the IR flux 
according to the SFRs in each of the two galaxies, as we have done 
for the LIRGs, the lower ~L IR in the Hwang et al. (2010) catalog 
yields two LIRGs. However, the slightly higher h IR in the Wang & 
Rowan-Robinson (2009) is sufficient to preserve the ULIRG clas- 
sification of the more highly star-forming component in the pair. 

Despite the lack of ULIRGs in our sample, it is worth consid- 
ering how the discussion in the previous section applies to ULIRGs, 
with log (L ;R /L ) > 12. To become a ULIRG, the ASFR has to be 
higher by a further order of magnitude than the diagonal dashed line 
in Figure[T3] which readily explains why ULIRGs are so rare in the 



local universe; unlike LIRGs, they are not 'naturally occurring' at 
any stellar mass. Even in mergers, the SFR boosts are rarely suffi- 
cient to produce a log (L IR /L Q ) > 12 in our SDSS sample. However, 
at high redshifts, the relationship between mass and SFR evolves, 
such that SFRs are higher at fixed stellar mass by a factor of 30 at 
z ~ 2 relative to z ~ (e.g. Daddi et al. 2007; Noeske et al. 2007; 
Wuyts et al. 201 1). This evolution means that a ASFR sufficient to 
produce a LIRG at z ~ (as in our SDSS sample) is sufficient to 
produce a ULIRG at z ~ 2. Indeed, we have (at most) 2 ULIRGs 
after SFRs are aportioned, ~ 0.5% of the sample, compared to 7% 
of pairs at z ~ 1 (Lin et al. (2007). 

There has recently been a shift in our understanding of the 
role of IR luminous galaxies and the implications of a classifica- 
tion scheme that has redshift independent h IR definitions. It has 
been noted in several works that high IR luminosities are the norm 
at z ~ 2 such that ULIRGs at low and high redshift are very dif- 
ferent beasts (e.g. Daddi et al. 2007, 2010; Elbaz et al. 2007, 2011; 
Sargent et al. 2012). Whilst a 'normal' M* ~ 10" M Q star-forming 
galaxy is a LIRG at z = 0, a similar mass galaxy will be a ULIRG 
at z = 2, and most galaxies with masses above 10 10 M Q will be 
LIRGs. Indeed, most ULIRGs are regular main sequence galaxies 
by z ~ 1 (Sargent et al. 2012). It is therefore not surprising that 
LIRGs dominate the cosmic star formation budget by z ~ 1 (Le 
Floc'h et al 2007; Magnelli et al 2009). Rather than using fixed 
limits in L /R , identifying the galaxies with the most enhanced SFRs 
at a given redshift may be best achieved by, as we have done here, 
quantifying offsets from the star-forming 'main sequence' of galax- 
ies at the same epoch (see also, e.g., Rodighiero et al. 201 1; Nordon 
et al. 2012; Sargent et al. 2012; Rovilos et al. 2012). Nonetheless, 
Kartaltepe et al. (2012) have argued that most z ~ 2 ULIRGS are 
still associated with mergers, despite apparently being in a steady, 
'maintenance' mode of star formation, rather than undergoing the 
intense starbursts asscoiated with the low redshift ULIRGs (Elbaz 
et al. 201 1). This could be explained if high redshift mergers tend to 
trigger much lower (if any) amounts of star formation than their low 
redshift counterparts. Such a scenario has recently been suggested 
by Xu et al. (2012) who do not find an enhancement in the SFR 
of mergers at z > 0.4 from Spitzer data. Determining the distribu- 
tion of SFRs amongst mergers relative to a well matched control 
sample is the natural next step in understanding the contribution of 
interactions to SFRs at high redshifts. 



5 SUMMARY 

From the SDSS DR7 we have compiled a sample of 9397 galaxies 
in close pairs with mass ratios 0.1 < Mi/M 2 < 10, projected sepa- 
rations r p si 80 hzl kpc and velocity differences AV < 300 km s -1 . 
The pairs sample has been matched to the SDSS-IRAS catalog of 
Hwang et al. (2010) in order to investigate the connection between 
mergers and high IR luminosities. Total IR luminosities are mea- 
sured for 400 of the pair galaxies, of which 172 are LIRGs and 1 
is a ULIRG. For comparison, a mass-, redshift- and environment- 
matched control sample of 4 galaxies per pair galaxy is also com- 
piled. Our main conclusions are as follows. 

(i) The LIRG fraction of galaxy pairs increases with decreas- 
ing projected separation relative to the control sample (Figure [4j. 
Such an increase is consistent with expectations that interactions 
are triggering star formation which, in turn, increase the IR lumi- 
nosity. Interestingly, we find an increased LIRG fraction out to at 
least 80 kpc (the limit of our sample), whereas many previous 
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Figure 13. ASFRs versus stellar mass for galaxies classified as LIRGs (large 
red circles) and non-LIRGs (small blue points). The diagonal dashed line 
shows the minimum ASFR required in order to produce a LIRG (from equa- 
tion[2). The top panel shows galaxies in our pairs sample; the control galax- 
ies are shown in the lower panel. 



LIRG studies have focused on morphological classifications and 
late stage mergers close to coalescence. 

(ii) Most of the LIRG pairs in our sample are biased to high 
mass ratios (Figure[5}> due to the higher SFRs that are intrinsically 
found at higher masses. The highest IR luminosities are therefore 
not necessarily found in galaxies with the highest SFR enhance- 
ments. 

(iii) The pairs with the highest L /R are generally those with the 
highest total mass, where the h IR reflects the total SFR in the IRAS 
beam. However, we do find that at fixed total stellar mass (M 2 +M 2 ), 
major mergers (0.33 < M1/M2 < 3) have, on average, slightly 
higher h IR than minor mergers (Figure |7}. The most likely candi- 
dates for (merger-induced) LIRGs are therefore close pairs where 
both galaxies are star-forming and have approximately equal, high 
stellar masses. 

(iv) LIRG pairs are found in relative low density environments 
compared with the full pairs sample. However, the environments 
of LIRGs are consistent with those of other star-forming (IRAS- 
detected) pairs. 

(v) The enhanced L ;R in the close pairs relative to the control is 
seen in both the star forming and composite (star forming+AGN) 
galaxies (Figure [8j. However, there is no statistical difference be- 
tween the h /R distributions of the pair AGN and control AGN. The 
enhanced h IR in close pairs is therefore likely to be driven by star 
formation rather than a higher incidence of AGN at close separa- 
tions. 



(vi) The AGN fraction increases with L ;R in both pairs and con- 
trol galaxies, although there is marginal evidence that the pair AGN 
fraction is lower than the control at fixed h IR (Figure[9j. 

(vii) The fraction of galaxies classified as close pairs increases 
with L; S above log (L; fi /L ) ~ 11 (Figure [TUt. Close pairs (r p < 30 

kpc, 0.33 < Mi /M 2 < 3) represent 7% of the LIRG population. 

(viii) LIRGs have higher SFRs and E(B - V) than non-LIRGs 
(Figures[TT1and[T2t. This is true for both pair and control galaxies. 
In fact, the SFRs and E(B - V) distributions are indistinguishable 
for LIRGs with/without a close companion. 

(ix) We define a SFR offset, ASFR, as the difference between 
a galaxy's SFR and the median of other star-forming galaxies at 
equivalent mass and redshift. The ASFR required to become a 
LIRG decreases as mass increases. A regular star-forming galaxy 
with mass above log (JVL/Mo) ~ 1 1 (at z ~ 0) is likely to have suffi- 
cient rates of star formation to be classified as a LIRG without any 
additional SFR enhancement. Lower mass galaxies require higher 
boosts; mergers can clearly contribute such boosts. 

The results presented here demonstrate the link between 
galaxy mergers and the high IR luminosities that define the LIRG 
population. However, we have also argued that interactions are just 
one (and, possibly, not the dominant) route to producing LIRGs. 
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